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RESEARCH SUMKARY 

Goldreich and Keeley continued their investigation of 
possible sources for the excitation of tne solar 5-minute 
oscillations. They also applied their linear non-adiabatic 
stability code to a preliminary study of the solar g-modes 
with periods near 160 minutes. In spite of a great deal 
of hard work, no definitive conclusions concerning the 
excitation of these modes were reached. Reports of the 
early stages of this work were presented at several conferences 
by Keeley and are included in the publications listed below. 

The excitation of the 5-minute oscillations by turbulent 
stresses in the convection zone remains a viable possibility. 
However, the power spectrum of the oscillations calculated 
by assuming that the turbulence is that given by mixing 
length theory does not agree very well with observation. 

In particular, the theoretical power spectrum declines too 
slowly with decreasing frequency and too rapidly with 
decreasing horizontal wave number to provide a good match 
to the observations. In spite of these difficulties, these 
calculations do give oscillation amplitudes which are of 
the correct order of magnitude. 

Theoretical calculations do not offer much support for 
the identification of the 160 minute global solar oscillation 
(reported by several independent observers) as a solar g-mode. 
The computations indicate that radiative dissipation in the 
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region below the convection zone exceeds the driving produced 
in the ionization zones by the tc-opacity mechanism. The 
large mass involved in these g-modes implies that even the 
substantial excitation they received from turbulent stresses 
produces a surface velocity amplitude which is much smaller 
than the observationally determined value. The most promising 
source of instability is nuclear driving, but this appears 
to be insufficient to overcome radiative damping. 

Our most recent work has been an attempt to reconcile 
mixing- length theory with the results of the calculations 
of linearly unstable normal modes . Here we have made a 
significant advance. Our calculations show that in a 
convective envelope prepared according to mixing length 
theory, the only linearly unstable modes are those which 
correspond to the turbulent eddies which are the basic element 
of the heuristic mixing length theory. We intend to write 
a paper entitled "In Defense of the Mixing Length Hypothesis." 

In addition to announcing the result just mentioned, we shall 
describe how the surface velocity field is related to that 
at depth in the convection zone. Ultimately, we hope our 
study will provide an explanation for the solar supergranulation. 
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j ,• Introduction 
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i . : 

| ■ Current. knowledge of tire interna! structure of the sun is almost entirely based 
! J upon theoretical calculations of solar models. Confidence in these models has 
; K;cn shaken by their failure to predict the solar neutrino fiux. The detcmiina- 
.* lion of the frequencies of the sun’s normal modes would provide us with sensitive 
♦ probes of the solar interior. If an appreciable number of low frequency modes 
**cre detected, the run of the basic variables p % p and T could be deduced with 
only a little help from theory. Even the internal rotation of the sun would be 
subject to observation study through the splitting of the modal frequencies it 
- produces. For the aforementioned reasons, recent claims of the detection of 
I. several long period solar oscillations have aroused considerable interest. We shall 
: review the main observational results and related theoretical investigations. 

* 

[ CLissijhation * 

' The normal modes of oscillation of the sun may be classified according to their 
** mgular eigenvalues / and m and the number of their radial nodes, w.f For each /, 

. dierc is u mode whose radial displacement and pressure perturbation have no 
: * udial nodes. This mode is designated the fundamental- or f-modc. Modes of higher 
frequency than the f-modc are called pressure- or p-modcs. The frequencies of 
die p-modcs increase with increasing u. For large n, the p-modcs resemble stand- 
{ . ing acoustic waves. The restoring force is picdominantly due to pressure, the 
velocity is approximately radial and is confined near the solar surface. For 1^=0, 
there is another sequence of modes whose frequencies arc lower than that of the 
f-mode. These are the gravity- or g-modes. The frequencies of the g-modes 
decrease with increasing n. For large /?, the restoring forces arc largely due to 

t For each l there are 2! + 1 nearly degenerate modes whose frequencies are 
* tightly split by the sun’s rotation 
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gravity which acts to reduce deviations of the surfaces of constant density from 
those of constant potential. The velocity is nearly horizontal and its amplitude 
peaks deep in the sun's interior. 1 

Observations of Excited Normal Modes 

Old results The power spectrum of the photospheric radial velocity field shows 
a broad peak (AP/P ~ 0.3) which is centered at a period of approximately 
minutes. The 5-minutc oscillations arc properly interpretated as the super- 
position ofmany, iion-radial p-niodcs. 2 ’ 3 ’ 4 Alt told there arc more than 10 6 
modes with periods between 4 and 6 minutes* For example, some combinations 
° f/ ( i nd ”’ denotcd b y (/.«)- which give p-modes with periods near 5-minutcs 

arc (0, 20), (10, 17), (I0 2 , 7) and (10 3 , 0). 

The solar granulation is the surface manifestation of turbulent convection 
Turbulent convection may be thought of as the superposition of many unstable 
iion-liucarly coupled, g-modcs. The non-linear couplings among these unstable ’ 

g-modcs arc so strong that the description of turbulent convection in terms or 
linear modes is generally not pursued. 


JVew results In recent years, several groups have reported the detection of 
additional modes of solar pulsation. It has been announced that the power 
spectrum of the temporal variation of the solar diameter exhibits several peaks 
at periods between 20 and 60 minutes. 6 Such periods fall in the range expected 
or p-modcs with small / and n. Indeed, some theoreticians have succeeded in 
matching the observed periods with those calculated front solar models 7 ’ 8 - 9 
If interpreted in terms of physical motions of the sun’s surface layers, the dia- 
l m f cr vanat *ons would imply surface velocities up to 30 m s*\ These inferred 
velocities are in conflict with upper limits of order 1 m s -1 established by direct 
observation of the photospheric velocity field. 10 It has been suggested that the 
apparent diameter oscillations are largely due to periodic variations in limb 
brightness rather than to surface displacements. 6 This suggestion is difficult 
to reconcile with the small observational limits placed on atmospheric tem- 
perature oscillations." Furthermore, we expect the sun’s surface layers to 
respond rather stiffly to driving at periods much longer than the acoustic cutoff ‘ 
peno o 1 80 s. This expectation is confirmed by detailed numerical calcula- 
tions of nonadiabatic p-modc eigenfunctions. These eigenfunctions show that 
variations in limb brightness make a negligible contribution to the apparent 
changes of solar diameter. 12 

Two groups claim to have detected an oscillation of the solar surface with a 
period of 2 hours 40 minutes and a velocity amplitude of 2 m s" 1 . ,3 ’ H The 
pulsation period is too long for any solar p-mode and can only be accounted for 
by a solar g-modc such as the g-mode with / = 2, n = 1 1 , 7 At first glance, the 
detection of this oscillation would appear to be totally convincing. Its discovery 
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»as announced simultaneously by two independent groups, and no reports of 
! conflicting observations have appeared in the astronomical literature. Ncverthe- 
| less, closer scrutiny raises some doubts. One of the groups based their claim of 
discovery on only 2 days of data. The period. 2 hours 40 minutes, is exactly 
|.«» ofa day. Finally, only a single mode was detected with a period longer than 
, I hour. Why is this mode so special since there arc many solar g-modcs with 
] periods longer than 1 hour? 

This brief description of the observational status of the detection of solar 
«>.rmal modes summarizes the published observations. It is not surprising cspcci- 
illy in a new and controversial subject such as ours, that rumors of additional 
i.vnilicani observations abound. However, in our opinion it is too early to draw 
any firm conclusions concerning the validity of the claimed detections of solar 
oscillations with periods hmger than 10 minutes. 

Since the observational situation is murky and seems unlikely to be 
dmlicd in the near future, we turn to thcoiy. that most unreliable ofastonem- 
imcal tools, to guide us. The theoretical goal is to estimate the energy in each of 
the solar normal modes. 


linear Stability Calculations 

Hoik in this area has largely been confined to studies of the p-modes. The prime 
Alliculty is the lack of any trustworthy thcoiy of time-dependent convection 
Nonadiabatic linear stability calculations which neglect the coupling between 
onvcction and pulsation predict instability for p-modcs with P > 5 minutes.*- 12 
•ic nving of the instability is due to the k - mechanism which operates in the H 
.omzation-H- opacity region just below the photosphere. The high frequency 
P-modes, P< 5 minutes, are stabilized by radiative damping in the solar 
ionosphere. The precise frequency that separates stable from unstable modes is 
wiy uncertain and depends sensitively on poorly determined details of the 
Jimosphenc structure. 

Unfortunately, the interaction between convection and pulsation does not 
ippear to be negligible. In fact, an investigation that includes the effects of time- 
dependent convection finds that all of the p-modcs are stable. 12 This investiga- 
ton uses a model of time-dependent convection based on mixing length theory, 
he modes arc stabilized by the mechanical damping due to turbulent viscosity 
m t ie convection zone. However, the calculated margin of stability is small and 
ik uncertainties in the theoretical model for turbulent viscosity arc considcr- 

Thus> 11 scerns wiscst to conclude that theoretical calculations are incapable 
,J determining the linear stability of the solar p-modes. 

One interesting corollary of the stability calculations is the determination of 
e quality factor, or Q, of the normal modes. The frequency width Aw of each 
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mode is related to the modal frequency cj and to Q by Aco/oo = 1/20. The ‘ 
frequency widths of the low order p-modcs arc much smaller than their frequency * 
splittings due to the solar rotation. For example, the radial f-mode has * ’ 

w= l-9 x l(T 3 s~ T ,@ = 2x 10 9 and thus Agj^4x 10" 13 s“ ! , while the rota- 
tion frequency of the solar surface ft 2.7 x 10‘ 6 s’ 1 . 

There is an indirect argument which suggests that the low order p- and g-modcs 
are not linearly unstable. It rests on the fact that the ‘observed’ amplitudes of 
these modes arc so small. The reported radial displacements A r/R^ s 10“ 5 . 

It is hard to imagine any non-linearities that might lead to the saturation of 
linearly unstable modes at these small amplitudes. 

If we accept the preceding argument and conclude that the low-order p- 
and g-modcs are linearly stable, we are still left with the problem of identifying 
their principal sources of excitation. Here we are guided by general thermodyna- 
mic principles. These tell us that sources of excitation are intimately linked to 
sources of damping. Since turbulent viscosity in the convection zone appears to 
be the principal damping mechanism for p-modcs with P>5 minutes, we are led 
to ask how convective turbulence might act as a source of excitation for these 
modes. 15 


Stochastic Fxcitation of the p-Modes by Turbulent Convection 

There arc two ways ol viewing the interaction between turbulent convection 
and a p-utodc. The turbulent convection may be described as the superposition 
ol many unstable g-modcs. The problem then reduces to one of calculating the ’ 
non-linear couplings between one p-modc .and two or more g-modes. The non- 
linear interactions among the g-modes broaden their frequencies so that resonance * 
relations required for energy sharing between the p-modc and two or more g- 
modes arc satisfied in many cases. Experience witli non-linear mode coupling 
in other circumstances leads us to expect the energy of each p-mude to be com- 
parable to the energy in those g-modcs which couple most strongly to t!u;t 
particular p-modc. This last qualification is necessary because the unstable g- 
modes in turbulent convection do not attain energy equipartition. Most likely, 
the g-modc energy spectrum has the KolmogorolT form. 1216 

The other way of estimating the excitation of the p-modcs by turbulent 
convection is more physical. Turbulent flows arc always accompanied by a 
fluctuating stress field which acts as a source of acoustic radiation. In an un- 
bounded medium, the dominant source term arises from the Reynolds stress 
and produces quadrupolc radiation. 17 The bandwidth of the emitted acoustic 
radiation spans the frequencies of the eddies in the turbulent flow. In the sun, 
acoustic disturbances (the p-modcs) arc restricted to a discrete set of frequencies. 
Turbulent convection in the sun can only generate radiation at the frequencies 
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,»f these discrete modes. The rate at which turbulent convection pumps energy 
into each p-mode depends upon the modal frequency and the spectrum of the 
turbulence. A plausible estimate of the turbulent spectrum may be obtained 
from a combination of mixing length theory and the Kolmogoroff scaling. ,J - 16 
Once the spectrum is chosen, it is a simple matter to calculate the rate of genera- 
tion of modal energy. The steady state value of the energy in each p-mode is 
jelcrinined by a balance between the power radiated at the modal frequency 
nul the mechanical dissipation of modal energy by turbulent viscosity. These 
interactions are such that the modal energy tends toward equipartition with the 
kinetic energy in a single turbulent eddy whose turnover time is comparable to 
the modal period. This simple picture is slightly complicated because the pro- 
perties of the turbulence vary with depth in the convection zone. It is easy to 
show that the strongest interactions between a p-modc and the turbulent con- 
nection take place at the level where the turnover time of the largest eddies 
fjuuincd to have linear dimensions of order the pressure scale height) is equal to 
ilie modal period. Detailed calculations based on a standard solar model yield 
values for the energies in the radial modes that range from 2 x 10 26 erg for the 
i-modc to 2 x 10 24 erg for the p 20 -mode. The absolute values of the predicted 
modal energies depend sensitively on the convective velocity. To a fair approxi- 
mation, the modal energies vary as the scv&nh power of the convective velocity 
(it fixed mixing length). Thus, until an accurate tF$oiy of turbulent convection is 
developed, it will not be possible to estimate tfee modal energies to better than a 
lew orders of magnitude. 

Once the energies in the p-modcs arc &;!atfrlcd, thciT surface velocity ampli- 
tudes may be predicted. The predicted velocities arc disappointingly small. 

Typical values for th radial f-through p 20 -mo?ies vary smoothly from 1(T 2 cm s'* 
to 6 x KT 1 cm $“*. The increase of the predicted surface velocity with n reflects 
the fact that as n increases the p-modcs become more strongly confined to the 
turfacc layers of the sun. 

There is no escaping the fact that the predicted surface velocities of the 
p-modcs arc much smaller than the values inferred from the observations of the 
tolar diameter. We see no hope of reconciling the theoretical and observational 
results. One or the other must be in error. 


the Five-minute Oscillations 

It has recently been established that the spectral power density of the 5-minutc 
oscillations is concentrated in bands in the plane. 18 * 19 Here k h Is the 
horizontal wave number and w is the angular frequency. The locations of these 
taids conform with theoretical predictions based on the assumption that the 
5-minute oscillations are the superposition of many excited, non-radial p-modes. 
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The precise locations of the bands are sensitive to the specific entropy in the 
convection zone. 20 The observations arc fit best by solar models which have 
convection zones of low specific entropy. Models of this type result if a high 
efficiency for sub-surface convection is assumed. These models have deep 
convection zones. 

The mechanism that determines the excitation of the 5-minutc oscillations is 
not firmly established. Several authors have proposed that (he p-modcs involved 
in (lie oscillations are linearly overstablc. 2 * 4,5 The driving of this overstability is 
attributed to the K-mcchauisni. Alternately, it is possible that the oscillations are 
stochastically excited by turbulent stresses in the convection zone. A crucial 
problem is the determination of the rate at which energy is shared among the 
p-modcs involved in the 5-minute oscillations. The transfer of energy among the 
p-modcs arises from their non-linear couplings. The most important couplings 
are within triplets of p-modcs. Once the rate of energy sharing among the p-modt 
is estimated (a problem we are working on), it can be compared to the rate at 
which the p-modcs share energy with the g-modcs. 

If the p-modcs couple more strongly to the g-modes than they do to each 
other, then the amplitude and spectrum of the 5-minutc oscillations would he 
due to the stochastic excitation of the p-modes by turbulent stresses in the con- 
vection zone. The high frequency cutoff of the spectrum could plausibly be 
attributed to the strong, radiative damping suffered by p-modes with P< 5 
minutes. Furthermore, we could reasonably conclude that the p-modes involved 
in the 5-minutes oscillations were linearly stable. In the absence of mutual 
energy sharing, we would not expect unstable p-modcs to saturate at the small 
amplitudes inferred from observation, nor would we expect the amplitudes to 
vary smoothly with / and n as observations suggest they do. 19 

On the other hand, if the coupling among the p-modcs is stronger than the 
coupling of the p-modcs to the g-modcs, the possibility that at least some of the 
p-modes are linearly unstable would have to be seriously entertained. If the 
p-modcs arc unstable and strongly coupled, they would form a system that i$ 
similar to that of the g-modes involved in the turbulent convection. 
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I. Introduction 


This paper presents a discussion of two of the important problems con- 
cerning solar oscillation modes. In- S II the possible effect of the mechanical 
boundary condition on the low frequency P modes is reviewed* and some limits 
to these effects are discussed in the framework of a simple analytical model. 

In- S III* a model for the excitation and damping of the 5 minute oscillations 
Is discussed and some results of calculations are presented. • S IV is a 
brief discussion of the global nature of the. oscillations. 

II. The Boundary Condition Problem for Long-Period P-modes. 

* ♦ < 

Hill's observations (1976) of long-period oscillations in the sun can 

be reconciled with low limits on the velocity amplitude if the amplitude of 
the temperature perturbation 6T/T is much larger than the displacement 
amplitude Sr/r at the atmospheric level of his observations. Observations of 
6T/T at lower levels of the atmosphere* in combination with Hill's inferred 
5T/T, lead to constraints on the gradient of 6T/T. The requirements are in- 
consistent with the results of conventional model calculation for the radial 
modes and low P modes* which have generally been assumed to be the most likely 
candidates. However* the model calculations are extremely sensitive to the 
boundary conditions* and it has been suggested that therein lies the resolution 
of the problem. Another possibility* which Hill has raised at this symposium, 
is that the oscillations may be g modes with spherical harmonic number L in the 
range 20-50. Preliminary calculations of these modes also appear to be incon- 
sistent with the observations. 

7 

a) Radial Modes 

Typical results of a conventional calculation for the fundamental 
radial mode are shown in figure 1. The fractional temperature pertur- 
bation is shown as a function of optical depth r* and the phase angle 
relative to the displacement amplitude is indicated at various points 


along the curve. The last point in the model cannot be shown on the 
optical depth scale , but the amplitude and phase are indicated in the 
upper left of the diagrata. The last plotted point is at an optical 
depth corresponding roughly to the temperature minimum. The normalization 

, iw 1 

of — is such that fir/r * 1 at the outer boundary of the model. Over 

the range shown, 6r/r varies by less than 1Z. For x > 1, the phase angle 

gradually approaches 180 s , and the amplitude approaches the adiabatic 

value by the point where the temperature of the model is about 20,000°. 

The temperature perturbation was calculated from 4M. 

©tC 

in which the flux and displacement, perturbations, and thus the real and 
Imaginary parts of th* eigenvalue, were calculated simultaneously. 

The mechanical bov&dary condition applied to the model was that 
the pressure perturbation vanish in the region outside the last boundary 
in the system of difference equations. It then follows that V* 6r in the 
last zone of the model is given by 

v-s? = %$rn( w * + ¥) Sr > 

in which p is the pressure, g is the . gravitational acceleration, and 
dm is the mass associated with the last boundary. The hydrostatic model 
was calculated with P ■ 0 in the exterior zone, and thus 

V-&? = #) Sr » 

where c is the soundspeed and H is the pressure scale height. 

A simple analytic model suffices to demonstrate the importance of 
the boundary condition. For simplicity consider a plane parallel . 
atmosphere with a constant pressure scale height. Then the displacement 
amplitude takes the form 

Sf — A &>ep(fc ir z') + B CXp(V” z ) j 

in which A and B are constants, and z is the height measured above some 


+ PVAT=-;l>V*F, 


arbitrary reference level. The wave numbers k 4 ’’ are given by 

= y "(/i O-ver)*). , 

2 2 

In which Y = 2H, and y* = '« /•( fig) • **»• waves are evanescent if k is 

real. The solutions fc* correspond to incoming and outgoir $ waves, respec- 

—3 

tlvely. For ui r> 2x10 , roughly the frequency of the- fundamental radial 

mode, y|Y«i if the scale height is much less than 10 10 cm. Thus, in 
the vicinity of the temperature minimum, where most model calculations are 
terminated, 

k + » H' r . k" fio’M/c* 

> 

The wave-number kw of the incoming evanescent wave is several orders of 
magnitude larger than k~ . The boundary condition used in the model 
calculation essentially amounts to allowing only the outgoing wave. (The 
4g/i term has been omitted because of the plane parallel assumption.) 

The point made by Hill et al. (1976) is that the adiabatic part of 
the temperature perturbation is given by 

= -crt-OfAP+BK’) 

if the point of interest is taken to be z - 0. Thus if the coefficients 
of the two waves are roughly equal in the displacement function, then the 
temperature perturbation is completely dominated by the incoming 
evanescent wave, and is very much greater than the displacement amplitude 
6r/r. Thus observations made in the vicinity of the temperature 
minimum might well show a large temperature amplitude. Lower in the 

atmosphere, the amplitude of the incoming wave has decayed by a factor of 

; 

the pressure ratio between the two points, since k”^ H , and the corres- 
ponding temperature perturbation will be much smaller. 

The question to be decided is whether in fact a significant amount of 
the incoming wave can be present. Such a wave can exist in principle 
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if there are any reflecting layers above the point in question. It can be 
demonstrated within the context of the simple analytical model above that 
for a single reflecting boundary with only an outgoing wave beyond it, 
|A/B|<1 at the boundary. This result is obvious from energy considerations 
if the region interior to the boundary allows the waves to propagate; it 
can also be shown simply for the evanescent case. 

In figure 2 the relationship between k and Y is plotted. Consider now 
a situation in which a boundary separates regions of different Y. Then if 
k* is the wave number of the outgoing wave in the exterior region, 

A C k -- k*) 

6 “ Ck*- M 

Both k“ and k* must lie on the lower curve, and k"*" is on the upper 
curve at the same Y value as k~. The curve 1/Y is halfway between the upper 
and lower branches. It is immediately obvious from the diagram that 
k - k‘ is bounded, whereas k - k* is not, so that if tha Y corresponding 
to k~ is not too close to the critical point, |A/bJ< 1 is easily satisfied. 
A little algebra shows that in fact it is satisfied for any choice of 
k* whatsoever, in the evanescent or propagating region, and for any pair k* 
in the evanescent region. 

In the case at hand, k*»k“, k*, and {A/8> f ^ k /k , I n 

Other' words, fc^A s k“B, which means that the two waves contribute roughly 
equally to temperature amplitude at the boundary. Therefore a single 
reflecting boundary near the temperature minimum will. never result in 
6T/T » dr/r. 




Unfortunately, It ie clear that a situation Involving more than one 
boundary cannot yield so clear-cut a result. Just as a coating on a lens 
can change the ratio of reflected to transmitted energy drastically, a 
similar situation occurs for evanescent waves.- Indeed, for suitable 
choice of scale heights in the exterior, middle, and interior regions one 
can arrange to have a purely incoming evanescent wave in the interior, a 
purely outgoing one in the exterior, and a mixture in the middle. Thus if one 
will admit the possibility of more than one reflecting boundary, the question 
Is still open. 

A calculation which goes fairly far towards closing it again is 

described below. An adiabatic eigenfunction calculation was performed on 

4 -3 

a solar model with a chromosphere extending out to T~ 2xi0 , P~ 10 

r* .• 

supplied by Roger Ulrich. The boundary condition applied was the same as 
before,. The result was a displacement eigenfunction in which 5r/r changed 
about 3 % from the temperature minimum to the outer point, and its derivative 
was essentially the same as that found in the previous model calculation. 

This demonstrates the important point that nothing drastic happens to the 
eigenfunction in the vicinity of the temperature minimum. That region does 
not act like a pair of reflecting boundaries, or at least, not a pair with 
serious consequences for the eigenfunction. The second point is that the 
region of the temperature minimum is seen to be separated from any other 
potentially-troublesome regions by at least six orders of magnitude in 
the pressure, over which any incoming evanescent wave would have decayed by 
roughly this* same factor. 

Thus it seems rather unlikely that the incoming evanescent wave can 
dominate the temperature perturbation amplitude at or below the temperature 


minimum. 
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The arguments above apply strictly to die adiabatic part of the temperature 
perturbation. In principle it is possible to do a full non-wliabatic calcula- 
tion of a model including a chromosphere, but in practice this can not be 
done without a good description of the heating and cooling processes, and 
their response to perturbations. An estimate of the possible effect of the 
radiative processes can be obtained in the following way. In a plane- 
parallel atmosphere the condition of no incoming radiation gives a relation 
between the energy density and the flux as t 0. If this relationship is 
perturbed about an equilibrium situation in which V*F (radiative) « 0, and 
applied at small optical depth to a model calculation, the dominant terms 
may be obtained from the condition 65 “■ a-THC » in which 

B is the Planck function and k is the opacity. This equation, along with the 

expression for $T/T given by the heat equation, leads to £?"= 

. v ■*" i.u>c.»t) > 

where i “ (-1/ and Cv is the specific heat in ergs/gm °K. Thus the total 

temperature perturbation is written as its adiabatic value divided by a 

factor which has magnitude > 1 in all cases* In the case of the solar model 

* . 

for which 6T/T was shown in figure 1, the factor in the denominator takes the 
value 1 - 4.341 , and the amplitude and phase it predicts for 6T/T are 
1.85 and -103°, as compared to 1.70 and -98° obtained by the model calculation 
in which some smaller terms were also included. 

In the presence of a non-radiative heating mechanism, the expression 
for 6T/T from the heat equation involves perturbations of the non-radiative 
contribution to the flux, and cannot be estimated without a more detailed 
discussion of such mechanisms. 

Besides the neglect of non-radiative heating and cooling mechanisms in 
the magnitude of 6T/T, another criticism could be directed at the treatment 
of radiative transport in optically— thin regions. Whether these holes in 
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current discussion are large enough to accommodate the temperature 
perturbation eigenfunction suggested by Hill's observations is unknown. 

It seems fairly likely that purely adiabatic contributions to 6T/T are 

Insufficient. 

b) Non-radial Modes 

The situation for non-radial modes is considerably more complex. The 
k(Y) relation for this case is Y^l ± (1 - 2y 2 Y + k^Y 2 )* 5 ) , in which 
k h « [L(L + lJ/R 2 ]* 5 is the horizontal wave number, and v 2m U Q 2 + kj, 2 /n o 2 . I^”\l- 1^ 
If 1 r “ 1 < u 2 /kt, < r,” 1 there is no propagation region for any value 

•* 1 O , n x 

of Y. The k(Y) relation is shown in figure 3. The reflection coefficient is 
much more complicated than in the radial case, and it is more convenient to 
deal directly with the coefficients of the Lagrangian pressure perturbation. 

If = A* exp(k + z) •+- B # e.xp(kTz) 

in the interior region, and 6P/P - C' exp (k"*Z) in the exterior, then 
continuity of 6P/P and the radial displacement at the boundary yield the 
result 



. 2 2 

Here cu^ 2 " c , and * refers to quantities outside the boundary, 

£ 

A simple case for which | A’ /B * ( >> 1 is obtained by choosing o) L ■ u. 

A* k* — k*’* 

The reflection coefficient reduces to •gr _ ^ 4 -vv) » It is clear from 

figure 3 that this ratio can be large. 
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For the case of p inodes with L < 20, model calculations including the 
chromosphere showed no sign of significant reflection in the region between 
the temperature minimum and the top point in the model. Higher L p-modes 
do not have frequencies in the range of interest.’ As in the radial mode 
case, the buffer provided by a factor of 10® in pressure between the last 
point In the model and the region of observation suggests that it is unlikely 
that there can be a strong reflected wave in the region of the temperature 
minimum. Of course, this conclusion has the same weaknesses as in the 
radial case. 

-3 

Gravity wave modes with L " 20, u ~ 2x10 have large amplitude even 
within 10% of the center of the sun, and complete solar models are required 
for a discussion of their nature. The chromosphere model supplied, by Ulrich, 
patched onto an interior model with a reasonably compatible equation of state, 
was used for an investigation of the g modes. 

The mode at (L ■ 20, tu » 2x10 ) had 16 nodes in 6r/r. (Modes with higher 
L at the same o* have more nodes and thus are increasingly difficult to cal- 
culate accurately by numerical techniques.) Both horizontal and vertical 
displacements varied quite slowly with radius above the surface of the 
convection zone and through the temperature minimum. Only in the immediate 

vicinity of the outer boundary point was there any significant contribution 

—3 *-3 

from the reflected wave. At (L * 30, u = 2.5x10 ) and (L ■ 50, u ■ 3x10 ) 

the behaviour of the eigenfunctions in the surface region was quite similar. 

The main difference came in the interior, where the amplitude was much 

higher relative to that at the surface; the result is that a kinetic energy 

orders of magnitude larger than is present in ail the 5 minute modes combined, 

would be necessary to produce a velocity of 1 cm/sec at the surface. 

These results suggest that the g modes are not significantly more favor- 
able than the p modes as an explanation of the observations. 


III. Excitation and Damping of the Five Minute Oscillations. 

The possibility that the observed oscillations in -the period range 
3 - 10 minutes could be excited by turbulence in the convection zone has 
been discussed by Goldreich and Keeley 1977. At the time of writing. 
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eigenfunctions for the non-radial inodes which dominate at these periods, were 
not available to the authors, so only a crude estimate could be made of the 

excitation to be expected. Extensive calculations of adiabatic eigenfunctions, 

/ ’ 

and the degree of excitation to be expected under various assumptions 

} 

about the solar convection zone have now been done; some of the initial 
results are presented below. 

One simple way of looking at the basic process is that the turbulent 
eddies generate acoustic energy over a broad range of frequencies determined 

I 

by the size and turnover time of the eddies, and their lifetimes as. 

Identifiable swirls. This acoustical energy constitutes an excitation of 
the normal modes of the sun as a whole. , Alternately, the turbulent 
convection may be regarded as an excitation of a large number of unstable 
gravity-wave modes, very strongly non-linear ly coupled to each other, and 
more weakly coupled to the P modes which dominate the 5 minute oscillations, 
and which are assumed to be accurately represented as linear oscillations. 

The latter approach, which treats the g and p modes on an equal footing, is 
more precise, but far less tractable than the former. Thus the first 

| 

picture was used as the basis for the calculations. 

Basically, what is done is to imagine that the properties of the turbulence j 

| 

are given, and not much influenced by the oscillations. The effect of the j 

j 

turbulence on the p modes is represented by inhomogeneous terms in the . 

I 

otherwise usual differential equations for the normal modes. The space and 
time behavior of the inhomogeneous terms is describable only in terms of 
the statistical properties of the turbulence. Conceptually, the problem is 
rather similar to that of Brownian motion of a particle in a harmonic oscillator 
potential. In the Brownian motion problem, the oscillator comes to energy 
equipartition with the molecules. The solar situation is more complicated 
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because it begins with a partial differential equation. However , in the 
limit where the turn-over time of the eddies is less than some numerical 
factor times the period of the oscillator (essentially the same as the 
assumption that the velocity of the Brownian particle changes negligibly 
on the timescale of the individual hits by molecules) the results obtained 
can be interpreted in terms of the normal mode coming to approximate energy 
equipartition with eddies having turnover-times approximately resonant with 
the normal mode. This bias towards one particular component of the 
turbulence presumably arises because of the assumption of a Kolmogoroff 
spectrum to describe the velocity at a given length scale; the kinetic 
energy resides mostly In the larger eddies. In the solar case, the dominant 
eddies In much of the convection zone (assumed to have size - scale height) 
have timescales much longer than 10 minutes. It is assumed that these- 
eddies are ineffective in exciting short-period modes. However, even at 
such places in the convection zone, there are smaller-scale eddies, with 
shorter timescales, which can provide effective excitation. Of course, 
because of their smaller size, they have less kinetic energy to share 
with the normal mode. 

The final amplitude of the oscillator depends on the balance between the 
energy added by the fluctuating forces, and the damping mechanisms. These 
are directly related, since the fluctuating forces which increase the kinetic 
energy of a normal mode from zero to the equilibrium value are also the 
origin of the frictional forces which reduce to the equilibrium value, the 
energy of a mode that initially is excited above that value. For this 
reason, the discussion necessarily turns to the dissipation of normal mode 
energy by turbulence. The excitation and dissipation can be treated 
together in the Brownian motion problem, but this has not been dme for the 


problem ae hand. ’ Instead, the effect of the turbulence has been modelled in 
terms of a turbulent viscosity coefficient. In keeping with the assumption 
made earlier with regard to excitation, it has been assumed that turbulent 
eddies with turn-over times greater than the oscillation period are ineffec- 
tive in providing a damping force. The main damping then comes from the 
largest turbulent eddies which satisfy the time-scale requirement; these are 
precisely the ones which are most effective in exciting the modes. A 
further complication, in ..the solar problem, which could have important 
consequences, is the presence of other damping and excitation mechanisms 
which also provide linear. damping in the equation of motion for the normal 
modes. 

The present calculations have been carried out using adiabatic eigen- 
functions, and the excitation energy has been calculated assuming turbulent 
damping only. The radiative contribution to the damping (or excitation) is 
currently being investigated. Implications of the radiative damping 
constants for the present calculations are discussed later. 

The eigenfunction and excitation energy are used to determine the 
surface velocity for each mode in the frequency range of interest. The 
number of normal modes involved is quite large. For example, between wave . 
numbers of .285 (Mm) * and .43 (Mm) ^ there are about 400,000 modes with 
periods in the range of 3 minutes to 10 minutes. It has been assumed that 
the velocities add incoherently, so the sum of the squares has been taken 
to determine the total. velocity. 

a) Results of the Calculations 

The rms velocity estimated earlier on the basis of knowledge of radial 
eigenfunctions only, was about a factor of 40 below the observed value of 
about 0.4 km/sec. With the non-radial mode calculation, the ordinary 


mixing-length theory with t/ H ** 1 gives 0.11 kra/sec. For t! H ■ 2, the 
result is 0.23 km/sec. Thus it still falls short by a factor of about 
4 in the energy required. 

The results are rather sensitive to the velocity predicted by the 
mixing-length theory. Let X be the scale and be the velocity of the 
eddies which make the dominant contribution to a given mode. Then the total 
energy of such an eddy is E * X 3 v^ 2 . The timescale ■ X/v^ * s to 
held constant as the mixing length velocity is varied , since the main 
contribution comes from resonant eddies. Assuming a Kolmogoroff spectrum, 
v^/v£ m (X/£) ^ 3 where t is the mixing length and is the convective 
velocity given by mixing-length theory. From these relationships it follows 
that E * V£^*^ if is held fixed. Thus a relatively small error in the 
convective velocity predicted by the mixing-length theory can make a 
large difference to the excitation energy, since the excitation energy is 
roughly proportional to t'u energy of the resonant eddies. 

Figure 4 shows the excitation energy of Mie lowest-frequency mode at 
a given horizontal wave number. The spherical harmonic index L = kR, where 
k is the horizontal wave number and R is the solar radius. The results are 
also shown if the mixing-length velocity is arbitrarily multiplied by a 
scale factor of 2. For the mixing-length unity model, the extra' fsr.tor of 
two leads to an rms surface velocity of .47 km/sec., and for the mixing- 
length two model, 1.36 km/sec. The range of periods corresponding to the 
wave numbers shown is 10 minutes down to about 4 minutes. The excitation 
energy as a function of frequency, for L ■ 300 (k 2 .43 (Mm) *) is shown in fig- 
ure 5. Each symbol represents a single normal mode, with the number of radial 
modes increasing with frequency. The large effect of an increase in the mixing 



length Is due in part to a higher density in the outer part of the convection 

rone, where the main contribution to these short period inodes arises. 

2 

The actual surface v is shown in figures 6 and 7. Note especially 

that although the excitation energy at a given wave-number is a decreasing 

function of frequency, the surface velocity Increases at first, because the 

kinetic energy of the higher modes is more favorably distributed. Finally, 

in figures 8 and 9, the observable energy is shown over a range of wave 

numbers and frequencies. Again, although the excitation energy of the high 

L modes is lower, much less kinetic energy is required for them to show a 

given surface Velocity. In all cases the surface point is at a low optical 
-3 -4 ' 

depth ”10 to 10 . If they were plotted at t * 1, then the high L 
modes, and the low L but high frequency modes would have somewhat lower 
velocities. Note that at L * 1500, only 2 different modes fall in the 
appropriate frequency range, whereas at L - 200 there are 11. This is more 
than compensated by the (2L + l)*fold degeneracy of the modes, so that more 
than half the total power is contributed by wave numbers above L ■ 1000. 

In general, the results are fairly encouraging for the model. Ulrich 
and Rhodes (1977) have suggested, on the basis of observations of the normal 
mode frequencies as a function of wave numbers, the mixing length parameter 
may be as high as 2. This would put the velocity predicted by the theory within 
a factor of two of the observations, without any tampering with the 
convective velocity. Since the theory predicts the energy as a function of 
frequency and horizontal wave number, it can be compared in great detail 
with observations such as Ulrich’s. Unfortunately, his complete set of data 
is not yet reduced, and the noise level in the part that is, is rather high. 
There is a possibility that the theory does not give a sharp enough increase 
in energy with increasing frequency, near 10 minutes period, - 
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There are several important loose ends* and limitations*' in the present 
calculations. Most of the excitation comes from the upper part of the 
convection zone where the timescales for the large eddies are not too 
long. This region in particular is likely to be sensitive to the convection 
theory. The velocities in the models used for the calculations presented 
here are surprisingly close to those in one of Ulrich's models* except for 
the outer scale height or so* in which they differed by as much as 50%. 

In all cases checked* less than about 25% of the excitation was produced in 
the discrepant region. 

The theory could easily be incorrect in the assumptions which were made 
concerning the spatial and temporal correlation functions* and the spectrum 
of £he turbulence. If the turbulent dissipation had been calculated directly 
by the same type of argument as the excitation rate* then there may have 
been more opportunity for errors to cancel. The turbulent viscosity coefficient 

was written simply as v^A, where X is the scale of the largest eddy which 

* 

has a turnover time < <•>***. Any extra numerical factor multiplying the 
viscosity* changes the excitation energy by the same factor* according to the 
present method of calculation. 

The calculations are incomplete because other damping mechanisms have 
not been included. The modes with periods shorter than about three minutes 
have been completely omitted because such modes are propagating rather than 
evanescent, in the vicinity of the temperature minimum. Thus they are 
expected to be damped quite strongly. The exact effect on the high-frequency 
end of the spectrum has not been’ calculated; this may be especially 
important at high wave-number, since much of the power is in those modes. 

Ando and Osaki (1975) calculated radiative damping rates, assuming that 
convection plays no role whatsoever. Unfortunately* the convective timescales 
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In the sun cover the range of periods of interest. Radial mode calculations 

Vith a simple treatment of time-dependent convection show that the convection 

may make the difference between a 'linearly unstable node, and a stable one. 

Calculations are currently underway to explore the situation for non-radial 

modes. Ando and Osaki find linear instability over a wide range of inodes 

involved in the 5 minute oscillations. This in itself is a little disturbing, 

unless some mechanism for limiting their amplitude to 6r/r ~ 10~* can be 
£ 

found. Unfortunately, the damping of the modes in this frequency range is 
very sensitive to the flux boundary condition (and possibly the treatment 
of radiative transport in the atmosphere). In the same approximation as 
discussed earlier, the perturbation of the flux divergence is given at the 
boundary point by + icJCyfJ m " AB ( I 3 ~ , 

where i ■ (-1)** For the frequency range of interest , the two terms on the 
left are comparable, and so the phase of 5(7*F/p) is sensitive to u>, and 
to the properties of the atmosphere. Although the radiative dissipation for 
T < 1 is not very large, the actual flux perturbation at t>> 1 is strongly 
influenced by the magnitude and phase determined by the boundary condition. 
Since almost all the radiative excitation takes, place in the layers at the 
top of the convection zone, the boundary condition has a very strong effect 
on the linear stability of the normal modes. 

It was noted that in the radial calculations of the mixing-length unity 
model, the turbulent damping rates were comparable to radiative excitation 
rates from a period of an hour down to 5 minutes. For the non-radial modes, 
this situation also holds, if the calculations of Ando and Osaki are reason- 
able estimates. In fact, there are modes which are slightly unstable in the 
linear sense, in spite of the turbulent viscosity. If the turbulent dam- 
ping is almost cancelled out by radiative excitation, then the steady-state 


* Hill has pointed out in this symposium that the .opacity perturbation 
SK/K is large enough for non-linear effects to be important when the 
temperature perturbation is summed over all the modes Involved in the 
five minute oscillations. This may be capable of increasing the rate 

at which energy is dispersed among the modes as a result of the non-linear 
coupling. 
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turbulent excitation energy of the mode will be much larger than the 
estimates given above. The turbulent damping rates of the mixing** 
length two model are not grossly different from those for mixing-length 
unity. However, those in which the convective velocity is artificially 
scaled up by a factor of 2 are roughly an order of magnitude higher. 

Direct: observations of the damping times of the modes would be very 
interesting, not only as a check on the radiative processes, but because 
of the relation of the turbulent damping to the excitation. Ulrich's 
observations show, for given wave-number, peaks at the expected normal 
mode frequencies. The widths of these peaks are uncertain because of 
the noise level, but it seems likely that it is due entirely to the spread 
of wavenumbers. Over the range from L ■ 10 to I> ■ 1500, the frequency 

h 

of the lowest mode at given I* varies very closely as L • Thus the expec- 
ted frequency spread is 6u>/u “ *s6k/k. The calculations indicate that 
the turbulent, damping increases with the wave number k. At L * 200 and 
300 the highest dissipation occurs in the modes with several radial nodes, 

but at L ** 1000 it decreases with increasing frequency. For the L = 1000 

-3 ‘ 

fundamental mode the "natural" width 6u»/tu • 3x10 if the convective 

velocity has its usual value, and is about 10 times larger if the velocity 
is scaled up by a factor of 2. The spatial resolution required to see the 
latter value is possibly within reach of ground based observations if 
guiding problems can be minimized and the seeing is very good. 

The problem of mode widths is discussed from a different viewpoint 
in the next section. 

IV Are the Oscillations Really Global? 

Do the solar oscillations know that they exist in a spherical body? 

If so, then the spectrum of horizontal wave numbers is discrete, and 


power In the k - u> plane should be found at distinct points rather than 
along continuous lines. 

The radial modes, the gravity wave modes, and the p modes with low 

horizontal wave number, all have appreciable amplitude deep down in the 

sun, and thus are expected to be "global'’ modes. However, the p modes 

with spherical harmonic index L > 200 are confined to a relatively thin 

surface layer comprising only 5% of the radius. The gross properties of 

these modes can be described quite adequately in a plane-parallel geometry, 

suc v . as Ulrich’s trapped wave model. Do these modes know, they are in a 

spherical shell rather than an infinite plane slab? Physically, the question 

is closely related to whether an inltially-locallized disturbance can 

propagate far enough in space and time to see the spherical geometry. 

If the damping time of a normal mode is T" 1 , it might be expected that if 

cy” 1 > circumference of the sun, the mode would be aware of the shape. 

However, the usual soundspeed may be far from the correct number to use • 

in this relation. In the present case, the dispersion relation is 

u 

known to have the behavior u> « k in the range of interest. Thus 
2Au>/u> ■ Ak/k - AL/L. If the mode width Au» is greater than <u/(2L), the 
dispersion relation allows a spread, in wave number corresponding to one 
unit in L. This is the damping required for the wave-number spectrum to 
be continuous. This gives as the condition for discrete modes, 

y < 3.15x10 * L ^ , for L - 300 the condition is 

Y < 1.82x10 This is less than a factor of two above the turbulent 
damping width for the 5 minute mode at I* - 300, when the usual convective 
velocity is used. Since the turbulent damping increases with L, this con- 
dition is likely to be violated somewhere in the range L ■ 300 to L ■ 1000. 
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Actually checking the discreteness of' the 5 minute inodes would be a 
difficult task. The individual modes sit in boxes with dimensions 
6k/k - 6u)/w s (21*) ** 1 in the k-u» plane. One or the other of these 

resolutions must be achieved if the question is to be decided. For 
periods in the range 3 minutes to 10 minutes, more than a day of observing 
time is required for this time resolution at L a 300 . The spatial resolu- 
tion is also difficult to reach. If the actual damping i3 higher than the 
estimate above, the transition from a discrete to a continuous wave number 
spectrum will occur at lower L where the time resolution problem is less 
severe. 

V. Summary of Conclusions 

1) It seems rather umlitely that low L p-modes can be responsible for 
the observed periods in the range 30-60 minutes. For g modes with L < 50 the 

current calculation^ also fail to show a significant contribution by incoming 

,» , > » 

waves to the adiabatic temperature perturbation near the temperature minimum. 
There are still some loopholes in' the theory, and it is not known whether they 
can make a significant difference. • 

2) The stochastic excitation model for the 5 minute oscillations, 
while very crude, is capable of giving a total excitation energy comparable 
to that which is observed, without much need to adjust possible parameters 
in the theory. The calculations at present do not give a prediction of 
the high-frequency end of the spectrum because of the uncertainty in radia- 
tive damping rates. 

3) " It seems fairly likely that the global oscillation model is not 
strictly applicable to the high wave-number modes involved in the 5 minute 


oscillations. 
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FIGURE CAPTION 1) The fractional temperature perturbation as a function 

qf optical depth. The fractional displacement pertur- 
bation is essentially unity over the entire range, and its 
phase angle is zero degrees. The phase of the temperature 
perturbation' relative to the displacement is Indicated 


FIGURE CAPTION 2) 


at several points along the curve. 

The wave-numbers k* as a function of scale height 

,.-l 


(Y s 2H) for radial modes. The curve Y lies half- 
way between the k + and k” curves. 

FIGURE CAPTION 3) The wave numbers k~ as a function of scale height, • 

for non-radial modes with horizontal wave number kjj. 

FIGURE CAPTION 4) The excitation energy of the lowest frequency p mode 

as a function of the spherical harmonic number L. 

Results are shown for mixing-length/ scale height » 1, 
and » 2, and for two choices of the convective velocity 
scale factor. A scale factor 1 corresponds to the 
ordinary mixing-length velocity. 

FIGURE CAPTION 5) The excitation energy as a function of frequency for 

the modes between about 3 and 10 minutes for 1* ■ 300. 

The plotted symbols denote discrete normal modes. 
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FIGURE CAPTIONS 6 and 7) The (surface velocity) predicted for the normal 

modes at l" 200 and b =■ 300. 
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FIGURE CAPTIONS 8 and 9) The (surface velocity) as a function of frequency 


and wave number, for two different mixing lengths 
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ABSTRACT 

The observed power spectrum of the solar five minute 
oscillations Is discussed from the viewpoint that the oscil- 
lations are excited by turbulent convection. The observations 
place significant constraints on the theory, and suggest 
constraints on the solar model structure. 


I. INTRODUCTION 




The solar five minute oscillations provide a sensitive probe of the 
properties of the solar atmosphere and upper convection zone, and a 
challenge to theories which attempt to explain the amplitude of the 
oscillations as a function of frequency and spatial wave number. This 
paper discusses some properties of the power spectrum of the oscillations, 

In the context of the specific excitation mechanism described by Goldreich 
and Keeley (1977). It also discusses the effect of the atmospheric 
temperature profile, and the mechanical boundary condition. In §111 the 
observed power spectrum as a function of aperture size is used to suggest 
constraints on the theoretical surface velocity of individual normal modes 
having periods near five minutes. It has already been noted by Ulrich and 
Rhodes (1977), that the frequency spectrum is better represented by solar 
models with mixing length equal to two or three pressure scale heights. 

It is shown in § IV that the steep low frequency side of the peak in the 
power spectrum is also more readily explained if the mixing length is 
greater than one scale height. The high frequency end of the power 
spectrum is also discussed. 

II. BRIEF DESCRIPTION OF OBSERVATIONAL DATA 
The main observational data to be considered here are the shape of 
the power spectrum, and the dependence of the power density at a given 
frequency, on the horizontal scale observed. Examples of the data available 
are given by Fossat and Ricord (1975), and Fossat, Grec, and Slaughter (1977). 
For observations through circular apertures, the power spectrum shows a 
rather sharp peak very near 5 minute period, and the position and shape are 
rather insensitive to aperture diameter over the range from 22" up to 
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several minutes of arc. One of the most striking features is the steep 
rise on the low frequency side of the peak; the power density increases 
by a factor of about 6.5 between about 7 minutes and 5 minutes. On the 
high frequency side, the drop-off is roughly half as steep. An additional 
piece of information is the fall-off in power density at the peak, as a 
function of aperture size, me data of Fossat and Ricord (1975) suggest 
a drop by a factor -2.5 in power density, when the aperture diameter goes 
from 22" to 60". 


III. THE RELATIVE CONTRIBUTION OF DIFFERENT SPHERICAL HARMONIC MODES 
a) Theoretical result of the averaging process. 

Consider observations made through a circular aperture sufficiently 
small that the region viewed on the solar surface subtends a small solid 
angle at the center of the sun. Then the sphericity of the surface can be 
neglected and the vertical component of velocity can be written in the 
form 

_ -iu>. t 

V(s. *, t) - 2 _ v Lnl Y Lra (8,»)e Lm , 

Lm 


in which e and <{> are spherical polar coordinate angles, Y, n) is a spherical 
harmonic, and tu is the oscillation frequency. v Lm is the velocity amplitude 
for a given Lm mode. The spatial averaging can be carried out most simply 
if the polar axis of the coordinates is chosen to be the line of sight. 

Let e c be the angular radius of the disc as seen from the center of the 


»■! 

Lr 

'i 


v 


Lm 



_[ 

Ml 


Lo 



t 


2n 

Ml 


680 


/ v u» da 

cos e o 

(ti 1 ) 7 **" Vl' < ms8 °>. 


4 


1 

*i 


sun. Then 


,-1 


where P^ Is an associated Legendre function. For L9 Q « 1, but 9*« 1 


V • I <^> ? 

L 


and for Le Q » 1 , but 0 Q « 1 , 



The theoretical calculations give mean square velocities for each 
normal mode; these are assumed to add Incoherently. Thus the formulae 
above yield 

V 2 \ L9 o« 1. and 

L 

V * *Z v l ^ cos ^ L 4 )e o ■ r] * L0 o >:> U 

To compare with observations made with a finite bandwidth Aw, the 
sum over L's is taken for all modes having frequencies within that band. 

If A u> is not too small, there will be one or more modes with no radial 
nodes, one or more with one radial node, one or more with two radial nodes, 
etc., contributing. Within each such group the L value will vary over a 
range depending on the bandwidth Aw, and the L values involved in different 
groups will be quite different except when L Itself is small or the bandwidth 
is wide. From the calculated relation between frequency and L for modes 
with a fixed number of radial nodes, 2 ~ . Thus the number of 

modes expected In Aw is —2 L ~ » provided this number is greater than unity; 
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* 

* 

If the bandwidth is very small some groups may not contribute at all. 

In this case, the bandwidth of the individual normal modes may be 

2 

important. If within each group it is assumed that V L is the same, the 
power can then be written as 

v 2 . L V 2 O < 2L >> Le o « 1 

v 2 ■ TT H * L 2 Sir Ttrp (2 L)cos 2 [(l4v I 1 ]. Le 0 »1 

L oo 

where the sum now is over the central L values of each group. The power 
per unit frequency interval follows directly from these formulae, 
b) Trial distributions of v^ 2 

It is instructive to plug in some trial distributions for v L 2 as a 
function of L, to see whether any clues to the actual distribution of 
mode energies in the sun can be obtained. The L values for the groups 
depend onto; for definiteness, a frequency u> = 2 x 10" 2 sec -1 was chosen. 
Approximate L values for groups near this frequency are as follows, for 
a model with mixing-length equal to one pressure scale height: 1004, 620, 

386, 270, 195, ; 50, 120, 98, 80, 68, 56, 50, 43, 36, 30, 25, 19. Relative 
power densities were calculated corresponding to the following three cases: 
Case 1 ) v L 2 = 1 for all L. Case 2) v L 2 = L. Case 3) v^-L" 1 . The 
calculations were simplified by omitting the cos 2 factor, and by using the 
asymptotic form for small L9 0 up to the point where it intersected the form 
for large Le 0 . This occurred at Le o « 2n"f . The results are shown for 
six different aperture sizes, in table 2. The most important point to note 
is that in all three cases, the ratio of power densities at 22" and 60" is 
greater than the observed ratio of 2.5 noted earlier. Although these 
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calculations are extremely crude, they suggest that there may be problems 
If v* Is constant or an increasing function of L. In fact, current 
calculations (Keeley 1977) suggest that v L a increases with L for periods 
longer than ~3.5 minutes, if turbulence provides the only damping mechanism. 
This suggests that it will be worthwhile to repeat the calculations 
presented in table 1, using an accurate representation of the Legendre 
function. Results similar to those of the crude calculation may present 
a severe challenge to the turbulent excitation theory, and a significant 
constraint on any theory which predicts amplitudes of individual normal 
modes. If, on the other hand, the observations at 22" suffer from seeing 
or guiding effects which reduce the power observed at high spatial wave 
number, then the results may be compatible. 

IV. THE SHAPE OF THE POWER SPECTRUM 

Since the observations suggest that the shape is relatively Independent 
of aperture size, it is convenient to discuss the power spectrum corresponding 
to fixed values of L. Results for L = 100, 200, and 300, all of which 
contribute to the power at periods as long as about 10 minutes, are 
considered in detail below. The discussion naturally divides into consider- 
ation of the energy to which an individual mode is excited, and the shape 
of the eigenfunction for the velocity amplitude. Of course, these are not 
totally independent, but this separation will be useful. 

The preliminary results reported by Keeley (1977) showed that for 
models with mixing length equal to one or two pressure scale heights, the 
surface (velocity) for 200 s Ls 600 had a peak near u> » 2 x 10“ (period 
~5.25 minutes). The peak was steepest on the low frequency side, in 
general agreement with observations. However, it was noted at that time 
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that the peak was not nearly sharp enough. These calculations include 
only turbulent dissipation in the calculation of the excitation energy, 
and the eigenfunctions were calculated for the adiabatic case. 

a) The shape of the eigenfunctions 

It is convenient to define an effective mass as the excitation energy 
required to produce a v 2 (averaged over the surface, and in time) of 
(1 cm/sec) 2 , at any particular depth in the atmosphere. This depends on 
the shape of the eigenfunction, but not on the actual excitation energy. 

The actual v 2 is obtained as the quotient of the excitation energy and the 
effective mass. 

On the low frequency side of the peak, the increase in v 2 is due to 
an initial decrease in effective mass as the number of radial nodes in the 
eigenfunction increases. Physically, this occurs because the kinetic 
energy of the higher modes is more concentrated in the surface region, 
where the density is lower, and less energy is required to produce a given 
velocity. In the models discussed by Keeley, the excitation energy decreased 
with to, but this effect was more than compensated by the decrease in effective 
mass, for to < 2 x 10” 2 . At higher frequencies (in a sequence with fixed L) 
the effective mass dropped off relatively slowly, and the net result was 
the high-frequency cut-off noted above. 

The problem of insufficient steepness at low frequency can be approached 
from two directions. The first is to construct models in which the fall-off 
of effective mass is more rapid in the frequency ranges = 1.5 x 10" 2 to 
2 x 10 -2 , and the second is to find models in which the excitation energy 
decreases more slowly withu), at least for periods greater than about 5 
minutes. The latter problem is discussed in b) below. The ratio of 
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effective masses over a given frequency range is conveniently expressed 
as - log (m(w 2 )/m(u>i))/log(”^-). In table 2 this ratio is shown for solar 
models with three different values of the mixing length. The eigenfrequencies 
are not the same In these three cases, but u>, « 1.5 x 10" and « 2 x 10~ 2 
for the functions chosen. It is clear from the table that the model with 
largest mixing length Is the most favorable, at all L values considered. 

If the excitation energies were roughly equal over this frequency range, 
the slope would be almost steep enough. Of course, the exact comparison 
with observations through a circular aperture requires that the results 
for various L values be combined as discussed in § 1 1 1 above. The results 
shown in table 2 are for models which have a fairly realistic atmosphere 
out to a temperature minimum of 4180°K. 

Some preliminary calculations of nonadiabatic eigenfunctions have 
also been done. Including turbulent viscosity in the equations of motion 
does not have a significant effect on the shape of the eigenfunction, as 
reflected in the effective masses for the low frequency, low L modes 
studied so far. On the other hand, a fully nonadiabatic treatment of 
the radiative dissipation has a significant effect, apparently because 
the dissipation is very strongly localized near the top of the convection 
zone . For the cases studied, the result is to steepen the decrease in 
effective masses with increasing frequency, and thus to steepen the 
low-frequency side of the peak of the power spectrum. The nonadiabatic 
calculations can't at present be done correctly with a realistic solar 
atmosphere, since the radiationflux is not given simply in terms of the 
temperature gradient, as in the approximation usually used for stellar 
interiors. Also, the effects of convection are not included. 
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1) Effect of surface temperature and the mechanical boundary condition. 

Models computed using the diffusion approximation all the way to the 
surface had surface temperatures of about 4880°. A more realistic value 
of T at the temperature minimum is about 4180° (Allen 1976). Adiabatic 
eigenfunctions and frequencies were computed for the diffusion models, and 
for models in which the empirical temperature profile was used at the 
surface. In addition, the models were computed with two different boundary 
conditions, one being that the lagrangian pressure perturbation vanish at 
the surface, and the other that an outgoing wave existed (evanescent or 
propagating) with radial wave number determined as if the surface had an 
isothermal region attached to it at the boundary point of the model. The 
effect of a lower temperature at the surface is to make the waves more 
evanescent, and is expected to be most important at high frequency. At 
periods & 5 minutes, it was found that with the 6^P = 0 boundary condition 
the low T model had lower effective masses (at optical depth 10" 3 ) and a 
slightly steeper decline in effective masses with increasing uj than the 
high surface temperature model. With the outgoing wave condition the 
situation was the opposite for both the magnitude of the effective mass, 
and the ratios of effective masses. In general, for either surface 
temperature, models with the outgoing wave condition had lower effective 
masses. The differences were largest at the highest frequency, where 
the waves were closest to being able to propagate. Significant changes 
in eigenfrequencies occurred only when the modes were close to propagating. 
Otherwise, the two atmosphere models and two boundary conditions gave 
nearly equal frequencies for the same physical oscillation modes. 
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b) The Excitation Energy 

In the theory described by Goldreich and Keeley (1977), the expression 
for the excitation energy of a normal mode is given in the form of a 
quotient of two numbers. The denominator is the damping rate of the 
normal mode* and the numerator is a double integral, over depth in the 
model, and over eddy sizes at a given depth. An Improved approximation to 
the inner integral has been used in the calculations described here. 

i) Low frequency behaviour 

In the discussion of the low frequency side of the peak in the power 
spectrum, it was noted that if the excitation energies over that frequency 
range were roughly equal, then the steep decline in effective mass for the 
high mixing-length models produced a slope much more in line with the 
observations. However, for all three mixing lengths tested, the excitation 
energy decreased significantly with increasing frequency. One way of 
equalizing the excitation energies for periods greater than about 5 minutes 
is to make all modes derive the main contribution to their excitation 
energy from a single set of eddies. The desired result will then be 
expected, provided the equipartition argument (Goldreich and Keeley 1977) 
is approximately valid. The result can be achieved, in fact, by a 
decrease in the correlation time for eddies in the outer part of the 
convection zone. If u> t < 1 at u> 4 2 x 10" 2 for the largest, most 
energetic eddies at a given depth, then all modes with w < 2 x 10” 2 will 
tend towards energy equipartition with these eddies. If u> r c > 1 for the 
largest eddies, then the modes will tend to equipartition with a smaller, 
less energetic eddy having a correlation time satisfying <u t 1 . In the 
present theory, the correlation time for the largest eddies is taken to be 
the mixing length divided by the convective velocity, and is scaled to 
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smaller eddies by assuming a Kolmogoroff spectrum. Then the Integral over 
eddy sizes depends on the correlation time assumed for the largest eddies 
In the form 



f & X 7 - 5 exp (-J-u,2t c 2 X2) 


J 

0 


This is a slow function of u> for tu t c « 1, but drops like w" 7 '® for 
u» t c » 1. The changeover between the two types of behaviour is rather 
gradual; thus a change In t which shifts the low frequency modes into 
the flat region results in a slower fall -off of the excitation energy >- 
and therefore of the power spectrum, at high frequency. 

The behaviour of excitation energy expected from the above discussion 
was verified by artificially increasing the convective velocity near the 
surface of the convection zone. A factor of less than two was sufficient 
to achieve the desired result. Thus it appears possible to reproduce, 
more or less, the steep rise in the power spectrum by a decrease in 
correlation time, in combination with a model with large mixing length. 

Of course, some of this gain is at the expense of the high-frequency 
fall -off. 

ii) High frequency behaviour 

The power spectrum as presently computed does not fall off very fast 
at high frequencies, at fixed L, except for L > 1000 or so (Keeley 1977). 
If such high L values do make substantial contribution to the power 
observed through typical apertures, they will improve the shape of the 
high frequency end substantially. However, it seems likely that radiative 
damping will play a significant role in decreasing the excitation energy 
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for high-frequency inodes. The only published nonadiabatic results on 
radiative damping are those of Ando and Osakl (1975, 1977). They find 
that the modes around 5 minutes are linearly unstable, but that modes 
near 3 minutes are stable. Strong damping sets In at somewhat lower 
frequency in their 1977 calculations, which include a chromosphere and 
corona. For the present models, the turbulent damping exceeds the 
radiative driving in all cases, but they are comparable for some modes. 

The result of a negative contribution to the damping would be an increase 
in the excitation energy of those modes; this could have a significant 
effect on the spectrum. For periods greater than about 7 minutes, the 
radiative growth or decay rate is relatively small compared to the 
turbulent decay rate. One important source of uncertainty is the 
calculation of the damping by the turbulent viscosity approximation. 

In addition, the calculations by Ando and Osaki do not include convection, 
and also do not calculate the radiative flux perturbation strictly 
correctly in the part of their model which employs an empirical T(t) 
relation. Some preliminary nonadiabatic calculations using models with 
the high- temperature boundary, but no convective perturbations, find 
stability for all modes checked in the 3 to 10 minute range, without the 
effect of turbulent viscosity. A further source of uncertainty is the 
convective velocity profile (and magnitude) in the outer part of the 
convection zone, since this region contributes strongly to both tht turbulent 
damping, and the total excitation. 
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V. SUMMARY OF RESULTS 


The crude calculation of the power density near five minute period, 
as a function of aperture size, Is not In agreement with observations; 

If a more accurate calculation establishes this discrepancy more firmly, 
then the observations may provide a powerful constraint on any theory of 
the excitation of the oscillations. In the present state of the theory, 
it seems possible to explain the steep low frequency side of the power 
spectrum peak. This requires, however, that the mass distribution of the 
sun be more like that of a model with a mixing length of three pressure 
scale heights, than one scale height. Important uncertainties in the 
damping due to radiative and convective energy transport preclude any 
strong statements about the high frequency end of the power spectrum; if 
there is linear driving comparable to the turbulent damping for periods 
near five minutes, this could have a significant effect on the sharpness 
and position of the peak. 

This work was supported in part by NASA grant NGL -05-002 -003. 
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TABLE 1 


Power density as a function of aperture size for three different assumptions about the velocities. 

Case 1: v^ 2 * 1 
Case 2: v L 2 = L 
Case 3: v^ 2 * L" 1 


Aperture 

diameter 

<2" 

5" 

10" 

22" 

60" 

120" 

e o 

< .00135 

.00325 

.0065 

.0143 

.0390 

.G78 

L c 

> 1004 

418 

209 

95 

35 

17 



Power density 

C 

c 


L < L c 

3.25 x 10° 

4.66 x 10° 

Case 1 

Power density 




L > L c 

0 

3.87 x 10 s 


Total 

3.25 x 10 6 

8.53 x 10 5 


Power density 

L < L c 2.69 x 10 9 1.84 x 10 8 



Power density 


Q 

Case 2 

L > L c 

0 

2.96 x 10° 


Total 

2.69 x 10 9 

4.80 x 10 8 


Power density 




^ < *"C 

6.52 x 10 J 

3.27 x 10 J 


Power der e Hy 



Case 3 

L > L c 

0 

5.32 x 10 2 


Total 

6.52 x 10 3 

3.80 x 10 3 


2.13 x 10 5 

4.38 x 10 4 

3.85 x 10 3 

~0 • 

1 .65 x 10 5 

6.83 x 10 4 

1 .33 x 10 4 

3.02 x 10 3 

3.78 x 10 5 

1.12 x 10 5 

1.72 x 10 4 

3.02 x 10 3 

2.97 x 10 7 

2.63 x 10 6 

1.01 x 10 5 

—0 

7.40 x 10 7 

1.39 x 10 7 

1.20 x 10 6 

-1.82 x 10 5 

1 .03 x 10 8 

7 

1 .65 x 10 

1.30 x 10 6 

-1.82 x 10 5 

1.95 x 10 3 

8.23 x 10 2 

1.51 x 10 2 

-0 ■ 

4.45 x 10 2 

4.68 x 10 2 

2.29 x 10 2 

8.7 x 10 1 

2.40 x 10 3 

1 .29 x 10 3 

3.80 x 10 2 

8.7 x 1C 1 




TABLE 2 

Logarithmic slope - ^^fog^y* ma - SS ^ between u» * 1.5 x 10" 2 and <u « 2 x 10" 2 , 
at three optical depths, for three values of mixing length/ pressure scale height. 



Mixing length 

1 

2 

3 


T - 10‘ 3 

2.596 

4.306 

4.944 

L ■ 100 

r » 10’ 2 

2.087 

3.929 

4.563 


t « 0.7 

1.2^3 

3.072 

3.688 


II 

O 

1 

U) 

2.846 

4.491 

5.060 

L - 200 

CNJ 

1 

o 

r— 

M 

2.426 

4.150 

4.719 


t » 0.7 

1.561 

3.378 

3.949 


T = 10‘ 3 

2.756 

4.746 

5.337 

L ■ 300 

T « 10" 2 

2.323 

4.419 

5.011 


T ■ 0.7 

1.517 

3.676 

4.271 
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EXCITATION OF SOLAR 6 HOOES WITH PERIOOS NEAR 160 MINUTES 


D. Keel ty 

Science Applications, Inc* 

Palo Alto, California 

ABSTRACT 

Solar g modes with t * 1 to 4 and periods near 160 minutes have been 
Investigated using a solar model with normal structure* Radiative dissipation in the 
region below the convection zone is much greater than the driving provided by nuclear 
reactions or the opacity mechanism. A crude treatment of convection suggests that it 
also is not an important source of driving. The damping due to turbulent viscosity 
is also small. Excitation of these modes by coupling to convective turbulence is 
substantial in terms of the rms energy of the modes, but the surface velocity is very 
small because of the large amount of mass Involved In the oscillation. 

1. INTRODUCTION 

The 2 h 40 m period first reported by Severny, Kotov and Tsap (1976) and by 
Brookes, Isaak and van der Raay (1976) has been supported by more recent observations 
(Scherrer et al. 1979). In this paper, properties of modes In this period range are 
examined in the context of a conventional solar model, and their stability 
investigated as described below. 

A model representing the sun was obtained by evolving a homogeneous 1 M 0 
model from the zero age main sequence until it resembled closely the present sun. 
The properties of the model were as follows: L * 3.83 x 10 33 erg sec"*, R * 6.9136 x 
10 10 cm, 2 * .02, Y(surface) * .244, X(center) * 0.39. The mixing length was 1.238 
pressure scale heights. The convection zone included about 1.2 x 10 31 grams and 
extended down to a temperature of about 1.5 x 10® K. 

The differential equations for adiabatic eigenfunctions were written in the 
form used by Osaki (1975), and solved by the method of Inverse iteration. The 
damping (or excitation) was then determined by perturbation theory, using the 
adiabatic eigenfunctions. The turbulent excitation energy and viscous damping were 
calculated as described by Goldreich and Keeley (1977), except for an Improvement in 
the Integral over the wave numbers in the turbulent spectrum. 

2. LINEAR STABILITY ANALYSIS 

The condition for linear instability is that a net positive amount of PdV 
work be done during one cycle of the oscillation. In calculating growth rates- 
according to perturbation theory, only the nonadiabatic part of the pressure 
perturbation can contribute to the work integral. If the time dependence is assumed 


r 


«K 

2<6 


to b« exp(lut), then the nonadlebetlc pressure perturbation is 


«P 


■ 

na 




( 1 ) 


In which 6 denotes a lagrangtan perturbation* and all other symbols have their usual 
meaning* For a given mass element, there Is a positive contribution to the driving 
If 



( 2 ) 


In which i Is the displacement. 

2*1* Radi a tive Contribution 

The equations used for the radiative flux ? r are 






(3) 


J - B « 


Tib 7>f v 


(4) 


In which J Is the mean Intensity, B Is the Planck function, and * Is the absorption 
coefficient. Since the perturbation of the scalar quantity ?• ? r Is required In 
equation (2), it Is convenient to write an equation for It directly: 



The second term on the left Is usually omitted In calculations of both the 
hydrostatic model and the damping. Its effect on hydrostatic models is normally 
small, but for some oscillation modes Its perturbations have a large effect on the 
damping. It is especially important at the cuter boundary of the solar convection 
zone. 

It is customary to do Eulerlan perturbations In nonradial calculations 
because they commute with spatial derivatives; this problem is more difficult to deal 
with in these circumstances than it Is for radial motions. However, In regions where 
* or other quantities vary very rapidly in space, the Eulerian perturbations are very 
much larger than Lagrangian perturbations. Because it is ultimately necessary to 
construct the Lagrangian perturbation for use in equation (2), the cancellation of 
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several significant digits, which could result when going from Eulerian to Lagrangii 

at the outer edge of the convection zone, could have a serious effect. Numeric* 

errors due to noise In the eigenfunctions could occur in this potentially importer 

driving region. Even if the eigenfunctions are not noisy, possible systemati 

effects due to the grid structure could Introduce significant errors. For thes 

reasons, Lagranglan perturbations were used throughout, and terms which wer 

cancelled In part or totally by similar terms from the convective flux perturbatic 

were handled explicitly. The final expression for the perturbation of v • = D i« 

EQUATION 6' r 

*0 - ij»CWD)ilc] sK = -5 -[b(b«B)5k]. k + jjs (iSD) ;I(J 

+ (t F r k) ;K " ( 5 J;K F r j);K * 5 J;K F rK;J < 6 ) 

in which 0 h (kp)" 1 . Index notation has been used because ordinary vector notatic 
is somewhat ambiguous for the last two terms. The semicolons denote covariar 
derivatives. 


2.2. Convective Contribution 

An equation for the time-dependent convective flux was written in the for 
given by Cox et al. (1966): 


< . 'c"’ • 

dt T 


(7) 


In which the timescale r is taken to be the mixing length divided by the loca 
convective velocity, and Is the convective flux as calculated from mlxir 
length theory. The perturbed form of the equation is 


4F c = T+Twt • (8) 


The Lagrangian perturbation of the instantaneous flux was taken to be 



in which a unit vector in the direction of -VP is taken as the direction of th 
convective flux, and is the magnitude of the flux as calculated in the usua 
way. Equation (9) expands out to the form 
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• f c l ” ■ !f C ,f • V ( H r ■ • «,>)»! • (10) 

In which H Is the pressure scale height, 

75 *& + ? 7 i * and *‘«r* v an + 5 i 7 i V ein 1 
The perturbation of the divergence of f c Is given by 

«(v.? e ) » 7.(6F C ) - v?:vF c , (11) 

in which the last term is the same form as In equation (6) for the radiative flux, 
and is partly cancelled by it. 

2.3. Nuclear Contribution 

The energy generation rate has been written in the usual form: 

e * c.p t ’ T v ■ (12) 

from which the perturbation in the adiabatic approximation is obtained: 

— - -V.?[n + (r 3 - l)v] . (13) 

Constant values for n and v were used throughout the energy-generating region. Since 
r 3 - 1 Is quite constant also, the nuclear contribution can easily be scaled to any 
desired values of the exponents n andv. For the results given In Table 1, the 
values n * 1, v ■ 15 were used. The high value for v was used because the nuclear 
reactions cannot maintain the equilibrium abundance of products In the PP chain for 
which the exponent v M or 5 is appropriate. Since nuclear driving did not appear 
to be an important effect, it was not considered necessary to treat this contribution 
In more detail. 

2.4. JM Condition for Local Driving 

With the approximations discussed above, the condition for driving of the 
Instability (equation 2) becomes 

7.1 1 + y + cv.t[n + (r 3 - l)vj| > 0 . (14) 

The nuclear reactions always contribute to Instability, as does the density 
perturbation In a region where v • F > 0. 
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Table 1 


TABLE 1, g Modes Near 2 h 4(f = 9600s 


£ * 

1 

2 

3 

4 

Period (s) 

10196 

9473 

9840 

9647 

Nodes 

6 

10 

15 

19 

Amplitude ratio 

12 

10 

22 

36 

Radiative damping 
rate (s' 1 ) 

1.40-13 

4.02-13 

9.16-13 

1.49-12 

Convective damping 
rate 

8.80-17 

-1.38-15 

-4.26-15 

-5.75-15 

Turbulent damping 
rate 

1.92-15 

1.50-14 

2.61-14 

2.85-14 

Nuclear 

excitation rate 

1.29-14 

1.33-14 

1.37-14 

1,38-14 

Radiative damping 
below convection 
zone 

99% 

98% 

96% 

96% 

Turbulent 
excitation 
energy (ergs) 

9.3+25 

1.6+26 

1.3+26 

7.9+25 

Energy for 1 cm/ sec 
at surface (ergs) 

1.1+33 

5.8+32 

8.7+32 

1.4+33 




3. RESULTS AND CONCLUSIONS 

Adiabatic eigenfunctions and frequencies were calculated for several modes at 
each value of £ In the range t « I to 4, and the work integral was calculated as 
described In 5 2 above. The results are summarized in Table 1 for the 4 modes with 
periods closest to 160 minutes. In addition, the radial and transverse components of 
X are shown In Figure 1 for the t * 1 mode. The 200 point grid resolved the 6 nodes 
of the £ ■ 1 mode quite well , but resol utlon was not as good for the 19 nodes at £ * 

4. However, it is unlikely that errors large enough to alter the conclusions are 
present. 

Except for the £ ■ 1 mode, the ratio of the maximum radial displacement to 
the surface displacement Increased strongly with t v as indicated In Table 1. 
Radiative damping was by far the dominant dissipative mechanism, and most of that 
damping occurred below the convection zone. Convective flux transport provided 
driving near the surface of the convecton zone in an amount comparable to the 
radiative driving there. The energy expected In these modes as a result of nonlinear 
coupling to convective turbulence was * 10 2 ® ergs, whereas the energy required for a 
surface velocity of 1 cm/sec is * 10 33 ergs. In any case, the turbulent excitation 
mechanism would not be a satisfactory one unless the observations eventually show 
that many modes are excited. 

The observations are not easily understood in terms of the conventional 
considerations discussed in this paper. 
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Figure 1» The radial and transverse components of the displacement vector for the i 
■ 1 mode. The eigenfunctions are normalized so that z r * 1 at the surface. The 
abscissa Is the grid point number. The lower boundary of the convection zone is at 
point 173. 




